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Abstract 

We develop a version of Haar and Holmgren methods which applies to 
discontinuous solutions of nonlinear hyperbolic systems and allows us to 
control the distance between two entropy solutions. The main difficulty 
is to cope with linear hyperbolic systems with discontinuous coefficients. 
Our main observation is that, while entropy solutions contain compres- 
sive shocks only, the averaged matrix associated with two such solutions 
has compressive or undercompressive shocks, but no rarefaction-shocks 
- which are recognized as a source for non-uniqueness and instability. 
Our Haar-Holmgren-type method rests on the geometry associated with 
the averaged matrix and takes into account adjoint problems and wave 
cancellations along generalized characteristics. It extends the method pro- 
posed earlier by LeFloch et al. for genuinely nonlinear systems. In the 
present paper, we cover solutions with small total variation and a class 
of systems with general flux that need not be genuinely nonlinear and 
includes for instance fluid dynamics equations. We prove that solutions 
generated by Glimm or front tracking schemes depend contiimously in the 

norm upon their initial data, by exhibiting an functional controling 
the distance between two solutions. 
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1 Introduction 



Our purpose is to develop a version of Haar and Holmgren methods which 
applies to discontinuous solutions of nonlinear hyperbolic systems and allows us 
to control the distance between two entropy solutions. The main difficulty 
is to cope with linear hyperbolic systems with discontinuous coefficients. Our 
main observation is that, while entropy solutions contain compressive shocks 
only, the averaged matrixv4(M, u') (see below) associated with two such solutions 
u, u' has compressive or undercompressive shocks, but no rarefaction-shocks - 
which are recognized as a source for non-uniqueness and instability, hence, the 
absence of rarefaction-shocks provides us with the key argument in order to 
establish the L}- stability of solutions. Our Haar-Holmgren-type method rests 
on the geometry associated with the averaged matrix and takes into account 
adjoint problems and wave cancellations along generalized characteristics. It 
extends the method proposed earlier by LeFloch and collaborators for the class 
of genuinely nonlinear systems [551 [5T] . 

In this paper, we present the framework in the context of the Haar method 
and cover solutions with small total variation and a class of systems with gen- 
eral flux that need not be genuinely nonlinear and includes for instance fluid 
dynamics equations. Our main contribution is a proof that solutions generated 
by Glimm or front tracking schemes depend continuously in the norm upon 
their initial data. This is achieved by exhibiting an functional controling the 
distance between two solutions. 

A completely different strategy is developed by Bressan and collaborators 
(see for instance [U [71 jTT] and the references therein) who rely on the concepts 
of generalized tangent vectors and shift differentials. This method leads to 
technical difficulties and requires either special properties of the hyperbolic sys- 
tems, strong regularity assumptions on the solutions, or sophisticated restarting 
procedure to keep the piecewise regularity of the (approximate) solutions. 

We discuss here the well-posedness theory for strictly hyperbolic systems of 
conservation laws 

dtu + d^]{u) = Q, u = -u(t,a;) e M^, (<,x) e M+ X M, (1.1) 

with general flux / : M^. In the present paper, we are motivated 

by a theorem established via independent methods in Bressan, Liu, and Yang 
[TUl [57] and Hu and LeFloch [25j for the class of genuinely nonlinear systems. 
Our main objective is to cover a class of hyperbolic systems (|1.1[) that need not 
be genuinely nonlinear and to establish that any solutions u, u' generated by 
the Glimm scheme or by front tracking satisfy the L} continuous dependence 
property 

|lw'(<)-«(<)|Ui(R) < ||7.'(0)-u(0)|Ui(K), teR+. (1.2) 

We consider solutions with small amplitude whose range is included in the 
open ball B^^ C K.^ centered at the origin and with radius So. In (|1.1[) . the 
flux / : B^,-, — !■ is a smooth mapping and, for all u, the Jacobian matrix 
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Df{u) admits N distinct eigenvalues Ai(m) < . . . < Xn{u) and, therefore, basis 
of left- and right-eigenvectors denoted by lj{u),rj{u) (1 < j < N). The latter 
are normalized so that right-eigenvectors have unit norm, and the basis of left- 
and right-eigenvectors are dual to each other. We refer to this as the "standard 
normalization" . 

It is well-known that solutions to (jl.ip are generally discontinuous and must 
be sought in the sense of distributions, typically within the class of functions 
with bounded variation. For the sake of uniqueness, the Cauchy problem asso- 
ciated with the system of conservation laws (jl.ip must be supplemented with 
an entropy condition, originally proposed by Lax [2 7) for genuinely nonlinear 
systems (satisfying, by definition, VA^ ■ Vi ^ 0) and later extended to general 
equations and systems by Oleinik, Wendroff , and Liu [39l |42l [34] . See the text- 
books [laillllllillST] for details. 

Pioneering work by Glimm [18] via the so-called random choice scheme estab- 
lished the existence of solutions to the Cauchy problem for genuinely nonlinear 
systems (|l.ll) and bounded variation (BV) initial data uq : M ^ B^^ : 

u{0,-) = uo. (1.3) 

In recent years, extensive research on nonconvex flux functions initiated [341 
[551 [231 [S [2 led to an existence theory for general systems [21 [SHI [3H]- In 
addition, Bianchini and Bressan 0] established the existence and continuous 
dependence of solutions to general systems (|l.ll) via the vanishing viscosity 
method; their technique does not cover solutions generated by Glimm's scheme 
or front tracking algorithms. Further important material on the existence theory 
can also be found in earlier works [40[ [5j [£[ [H [9l [211 [SI] on genuinely nonlinear 
systems. 

To establish the continuous dependence property (|1.2p we follow [25l [31] and 
introduce a matrix- valued map A =A{u, u') such that 

A{uy){u' ~u)^f{u')~f{u), ^^^^ 
A{u,u') = A{u' ,u), u,u'eB5o. 

We refer to A{u^ u') as the averaged matrix of u, u' and, for instance, we could 
choose ^ 

A{u,u')::^ [ Df{{l-e)u + eu')d9, (1.5) 
Jo 

although a different choice may be preferable for certain systems (jl.ll) . for in- 
stance to handle solutions with large amplitude. Since / is strictly hyperbolic 
we can assume (for a smaller value of Sq, if necessary) that^(M, u') admits N 
real and distinct eigenvalues denoted by 

Xi{u, u') < . . . < Xn{u, u'), 

together with left- and right-eigenvectors lj{u, u'),rj{u, u'). 



3 



If M = u{t,x) and u' = u'{t,x) are solutions to the nonlinear system (jl.ip . 
then clearly tp := u — u' is a solution to the linear hyperbolic system 

9tV + a4AV)=0, (1.6) 

where A := A{u,u'). Consequently, in order to establish ()1.2p it is sufRcicnt to 
derive the stability property 

\m)\\LHR)<um\LHM), teR+, (1.7) 

for a sufficiently large class (covering the situation of interest) of matrix-valued 
mappings A = A(f, x) and solutions ip = tp{t^ x) (which need not be of the form 
A =A{u^ u') and = u — u'). 

One should realize that the stability estimate (|1.7p does not follow from the 
standard theory of linear hyperbolic equations, since (|1.6p has discontinuous 
coefficients. In fact, we should not expect the estimate (jl.Tp to be valid for 
arbitrary systems (|1.6p : see Example 12.21 One of our objectives is precisely 
to identify a sufficiently large class of systems (|1.7p covering the situation of 
interest and, especially, to identify properties that the averaged matrix A(u,it') 
must satisfy. 

In most of the present paper, we assume that ip and A are of the form 
Ip = u — u' and A := A(u,u') and postpone the discussion of more general 
solutions and matrices. In practice, we will obtain (|1.2p (or (|1.7p ') only in the 
limit of a sequence of approximate solutions Uh, u'f^- 

\\ui,{t)-u'^{t)\\Li(t,)<\\uh{0)-u',M\\LHR)+o{l), t>0, (1.8) 

where o(l) tends to zero with the discretization parameter h. 

The approach above was used by Haar [22 to establish uniqueness results 
for smooth solutions of partial differential equations. It is also formally equiv- 
alent (by duality) to Holmgren method. The method was applied to nonlinear 
hyperbolic equations with special structure by several authors, beginning with 
Oleinik [39 for scalar equations. (Cf. [Ml [33] and the references therein.) Us- 
ing this strategy to handle discontinuous solutions of general nonlinear systems 
(|l.ip was first proposed by LeFloch and successfully led to a proof of ()1.2p for 
genuinely nonlinear systems ^25] . Independently, a method based on an explicit 
weighted distance was introduced by Liu and Yang [35] for scalar conserva- 
tion laws with convex flux and later extended to genuinely nonlinear systems 
by Bressan, Liu, and Yang [lOl [37] . 

The content of the present paper is as follows. 

1. Recall that a rarefaction-shock in the i-family, by definition, has all i- 
characteristics moving away from it. We prove that rarefaction-shocks are 
the only source of non-uniqueness of solutions to (|1.6p and provide an 
estimate for piecewise constant solutions of ()1.6p : 
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where the error term I'tja (t) involves only rarefaction shocks in the matrix 
A. 

2. Generalizing a similar observation made earlier for genuinely nonlinear 
systems we prove that the averaged matrix u') associated with 
two entropy solutions u, u' does not admit rarefaction-shocks. More pre- 
cisely, this property holds for "sufficiently robust" shocks, in a sense de- 
fined below (Section[3]). Together with our stability estimate, this gives 
a strong indication that (|1.7[) should hold. 

3. Next, we actually construct an functional satisfying the condition re- 
quired for stability. Our new functional can be regarded as a generalization 
of functionals proposed independently by Liu and Yang (via an explicit 
formula) and LeFloch (via a constructive method) for genuinely nonlinear 
systems. The main novelty is that not only the sign of (the characteristic 
components of) u' — u must be taken into account but also the monotonic- 
ity of the waves (the sign of u+ — m_ or tt^ — u'_ at jumps) and the jump 
of the averaged speed Ai(u+, u'_^) — Ai(w_, u'_). Indeed, the eigenvalues of 
the averaged matrix A(w, u') and their jumps across shocks are found to 
play a central role (instead of the usual wave strengths) . 

4. Furthermore, our construction of the weight takes into account wave can- 
cellations, by solving a backward evolution problem for an adjoint equation 
associated with (|1.6|) . Wave partitions (Liu [35]) are essential to exhibit 
the cancellation effect. Our present approach is in fact a generalization of 
the construction of the weight proposed in [25l [31] for of genuinely non- 
linear systems. (Namely, trains of oscillating waves with alternating signs 
and equal strength already arose in ^5] I31j.) 

5. We conclude with the continuous dependence property (|1.2p for a large 
class of hyperbolic systems, including the equations arising in isentropic 
fluid dynamics. 

Our method should be useful in other applications when wave cancellations 
allow to derive sharper bounds on solutions. For the Lagrangian and Eulerian 
compressible fluid equations, our method applies to solutions with large total 
variation. We refer to a follow-up paper for further developments. 

2 Shock waves and characteristic components 

2.1 Compressive, undercompressive, and rarefaction shocks 

We will rely on the following terminology and notation. In this section and 
in Section 2] the maps ijj = ^p{t, x) and A = A(t, x) are piecewise constant and 
admit finitely many polygonal lines of discontinuity, which can intersect two at 
a time and at finitely many times only; we refer to these lines as shock waves. 
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Consider, for instance, the vector-valued map A, together with the associated 
decomposition 

K+ X M = C^ UJ^UJ^, 

where is the set of points where A is locally constant, the set of points 

where A has a shock, and the set of points where two shock waves meet. 
Also, J^{t) will denote the slice of consisting of jump points with time 
component t. At each jump point {t,x) e J^, denote the left- and right-hand 
limits of A by A* = A*(t,a;) and the shock speed by A = A {t,x). 

We assume that, at each {t,x) € C^, the matrix A{t,x) admits N real and 
distinct eigenvalues 

Af-(i,a;) < ... < A^(i,a;) 

and, therefore, left- and right-eigenvectors l^{t,x),r^{t,x) (1 < J < A^), nor- 
malized in the standard way (see introduction). We always assume that the 
basis if' and rj^ remain sufficiently close to fixed basis of constant vectors. At 
each {t,x) e J-^, the notation Ij^ = Zj^(La;) and rjj_ = rjj_(t, a;) for the 
corresponding left- and right-hand limits will be used. 

We also assume that the eigenvalues are uniformly separated, in the sense 
that there exist fixed and disjoint intervals [A™™, A™*"] (j = 1, . . . , A'') so that 
at each {t, x) e C'^ 

A""" < Xfitx) < Xf^"", 
while for each {t,x) G there exists an index i such that 

Xf" <X^{t,x) < A^^. 

This allows us to decompose the set of shocks in A, in form = J^^L). . .Uj^, 

where J'/^ is the set of all i-shocks. 

Our results will be uniform for all matrices satisfying the uniform total vari- 
ation bound 

TV{A{t)) <C, t>0, (2.1) 

for some fixed constant C > 0. A matrix-valued field A = A(t, x) fulfilling 
all of the conditions above will be called uniformly hyperbolic with bounded 
variation, and we will be interested in estimates that are valid for the whole 
class of such matrix-fields. 

Definition 2.1. An i-shock wave in is called 

1. a compressive (or Lax) shock if Xfi_ > A"*^ > A-^, 

2. a slow undercompressive shock i/min(A^, A^) > A , 

3. a fast undercompressive shock i/max(A^, A^) < A"^, 
4- or a rarefaction shock if Xfi_ < )^ < A^. 
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Based on this definition, a shock which, for instance, is characteristic on its 
left-hand side and compressive on its right-hand side can be regarded as either 
a compressive or a fast undercompressive shock. For convenience in certain 
formulas, we may also allow large inequalities in the definition of rarefaction- 
shocks. 

The notation S^, J--^, and TZ-^ will be used for the sets of all com- 
pressive, slow undercompressive, fast undercompressive, and rarefaction shocks, 
respectively, and Cf-, Sf^, Tf^, and TZf- will denote the subsets associated with 
the ^-characteristic family. Hence, 

Example 2.2. It is not difficult to see that when N ~ 1 and A — a : K+ x M 

R is chosen to be 

then the Cauchy problem associated with the scalar equation 

dtip + d^ia^;) =0, = i^it, x) £ R (2.2) 

admits infinitely many solutions. Therefore, if the speed coefficient el admits 
rarefaction- shocks then the ( uniqueness property and, a fortiori, the ) continuous 
dependence property (|1.7[) can not hold for solutions of (j2.2p . 



To establish (|1.7p one possible strategy is to determine a weighted norm 
that is decreasing in time for all solutions of (jl.6p . In view of Example [2T2l above, 
we will need to restrict attention to matrix A with no rarefaction-shocks. In- 
troduce the characteristic components a — (ai, . . . , un) of a solution ■(/;, defined 

by 

=: ^aj(i,a;)rf {t, x) e n C''' . (2.3) 

j 

All the summations are over {!,..., N}, except if indicated otherwise. Then, 
consider the weighted norms 



l'/'(Ollw(t) := / V |aj(t,a;)| Wj(i,a;)da;, 



where the piecewise constant weight (wi, . . . ,WAr) satisfies 

w'"'" < wj(i,x) < w'""" (2.4) 

for all j and {t,x). Clearly, for fixed constants u;™!" ^^max ^ these norms 
are equivalent to the standard norm. 

Following [25l [31] we require that w be a solution of the adjoint system 
associated with (|1.6p 

dtMv + Ad.j:W = g, (2.5) 
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where the source-term g may consist of measure-terms and will be specified later 
when we will impose additional constraints on the jumps of the weight. (Strictly 
speaking, at a shock of A a definition of the nonconservative product A dx^ff is 
necessary, and indeed a suitable definition can be given along the same lines as 
in [551 1301 HSl 1311 ■ In the present paper we bypass the difliiculty and need not 
address this issue directly.) 

Precisely, we assume that every shock in w coincides with either a shock in 
A or a characteristic straightline, so that for all but finitely many (interaction) 
points (i, x) £ we have either 

1. {t,x) e J^, or 

2. there exists an index i such that A (t, x) = X'^{t,x) and, for all j ^ i, 

^) = wj(t, x). 

In addition, we could also allow the components Wj to exhibit a decreasing jump 
in time, that is, 

■Wj(t+,x) < Wj{t-,x), {t,x) e M+ X M. (2.6) 



We will say in short that w formally solves the adjoint system (|2.5p and, later. 



The following two lemmas were established in ^25. ,31^ 



will also impose further constraints on the jumps — w^- 



Lemma 2.3 (Time-derivative of the weighted norm). Suppose that the weight 
w formally solves the adjoint system. Then, for any piecewise constant ip the 
weighted norm satisfies 

j^\\m\Ut)<J2 E Piit,x)^j{t,x)+p+it,x)^+it,x) (2.7) 

for all but finitely many t, where 



p-{t,x) := (A {t,x)-\f_{t,x)) \aj{t,x)\, 
p+{t,x) (Af+(<,a;)-A^(t,x)) \a+{t,x)\. 
are referred to as the characteristic flux. 



(2.8) 



Lemma 2.4 (Signs of the characteristic flux). Across each i-shock the compo- 
nents satisfy for all j ^ i 



and for j = i 



±/3f<0, j<i, 
±/3f>0, j>z. 



< 0, C^, Pt > 0, TZf, 



(2.9) 



(2.10) 
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Clearly, the signs (and amplitude) of the terms A {t, x) — Xfj. {t, x) are es- 
sential in the evolution of the weighted norm. For j i one of these two terms 
always has a favorable sign. For j = i the signs depend upon the nature of the 
shock, and it is only in the case of a rarefaction shock that none of these two 
terms have a favorable sign. 



2.2 Properties of characteristic flux 

Our estimates below depend upon the strength of i-shocks in the matrix A = A(t, x), 
measures with the total strength |A+(i,x) — A~{t,x)\, or in term of the i- 
eigenvector \rf^{t,x) — rft{t,x)\. It is convenient to introduce also the charac- 
teristic flux with the modulus suppressed in (|2.; 



Lemma 2.5 (Jump relation for characteristic flux across an i-shock). Let ip = 

tp{t,x) he a solution of the system (jl.6p . Then, at every shock point {t,x) the 
characteristic flux 

,f := (a" - A- ) of 

satisfy 

lt-lJ+T.^fArt-rt)lk (2.11) 



k 



and 



^3 



. '\ ■ E (rt rt) + (7.7 - 7.+) r^)- (2.12) 

Lemma 12.51 implies the expansion 
lt^lJ+0{\rt-rt\)\ir\ + 0{\A+^A-\)J2W-\, l<j<N, (2.13) 

and for j i and j = i, respectively, 

7t-l- + 0{l)Y,\j+-^^\+Oi\rf^^-rt\)\ir\+Oi\A+-A-\) ^ \a^\, 

k^j k=ii,j 

7+ = 7r + 0(1) E |7,+ -lk\+ 0(|A+ - A-|) ^ \a^\. 

(2.14) 

Observe that (|2.14p contain no error term in the indices j and z, respectively. 
For instance, the second identity shows that 



7+ = 7r + 0(l)E 



k I' 



± 



and, therefore, if a ■ = oi^ = for all j ^ i then = •y- 
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Proof. The Rankine-Hugoniot relation associated with (|1.6p reads 



thus 



k 



(2.15) 



Muhiplying this identity by each eigenvector and using the normahzation 
lf+ = if j 7^ fc and / rf_^ = 1 , we arrive at (|2.1ip . 

For iUustration, suppose for instance that, for some fixed j, aU of the com- 
ponents 7^ for k ^ j vanish. Then, the identity (|2.15p imphes that 

Hence, either 7" = 7^" = 0, or else r-^_ = rj^ (since both are unit vectors and 
are close to each other) and, in both cases. 

We observe that there should be no error term (for instance in (|2.13p ) when the 
components 7^ vanish for all k but one single index. 

We can reformulate the estimate (|2.1ip by observing that for all j 



+ A — A \ ^ — A — A 

^] - 7j ^3- =2^1k - Ik ^k- 



= E 7.7 {rt ^rt) + ilk - 7,+ ) ^■ 



(2.16) 



Since the 's are unit vectors 

7/=77rf+.rf_+rf+.r!, 
77=7+^--^-^- -17, 

this implies 

(1 - (rjV • rfj^) 7; = - (rf+ • rf_)rf_) ■ n, 
and, since rj^ • rj^ = 1 + o(l), 

Combining this estimate with the second identity in (|2.16p we obtain (|2.12p . □ 
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We fix a (sufficiently small) constant ki > and introduce the following 
terminology. 

Definition 2.6. Consider a solution tp — ijj(t, x) of (jl.6p together with its 
characteristic flux (3 j defined in (|2.8p . Then, at any i-shock the j -characteristic 
flux {I < j < N) is said to be dominant if 

^1 |/3-| > \rt rt\ |/3-| + |A+ - A^j (2.17) 

The notion of dominant flux is quite natural in view of the expansion (|2.13p . 
Indeed, it follows that for all dominant characteristic flux 

sgn(7/) = sgn(7^^), j dominant, (2.18) 

which leads us to the following conclusion. 

Lemma 2.7 (Signs of characteristic components). For all j ^ i 

sgn{a^) — sgn{aj), j dominant, (2-19) 

while for j — i 

/ j-N I sgn(ay), CiUTZi and i dominant, 
sgn{aj) — < (2.20) 
\~sgn[a^ ), SiUJ-i and i dominant. 

Remark 2.8. 1. From the expansion (j2.13p and the inequality (|2.17p . an anal- 
ogous inequality follows with (ij replaced by (3^, i.e. 

2m 1/3+ 1 > |rt - rt | |/3r I + I A+ - A" | ^ \f3^\. 

k 

Conversely, if the latter inequality holds, then f3~ satisfies (j2.17p with ki re- 
placed by 3ki . 

2. At this stage of the discussion, in view of Lemma \2.4\ it should be clear 
that for the right-hand side of (j2.7p to be non-negative the weight should jump 
up or down at i- shocks, according to 

'>0, j<z, 

< 0, j = i and S^, 

> 0, i ~ i and Ti, 

.<0, j>i. 



w+ - w^. 



3 The averaged matrix admits no rarefaction- 
shocks 

3.1 Scalar and fluid dynamics equations 

Our aim in the present section is to derive properties satisfied by the averaged 
matrix A(u, u') (see (|1.4p and, for instance, ()1.5|) ) associated with two solutions 
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u,u' of Recall that a shock is entropy admissible if it satisfies the stan- 

dard entropy condition [Ml [13 [42l [34] ; see (O and ([3:26|) below. Our main 
observation is that entropy admissible shocks in weak solutions u, u' can not 
be rarefaction-shocks for the matrix-field u'). Throughout this section, we 
consider weak solutions with bounded variation, which are known to admit, in 
a measure-theoretic sense, shock curves and left- and right-hand traces. 



Precisely, we consider maps u in the space (R+ xl 



)nL'= 



pN 



of (bounded) functions with bounded variation. To such u we can associate a 
decomposition of x M (e.g. Federer [T7] ) 



(3.1) 



where C" is the set of points of approximate continuity of u, J'^ the set of 
approximate jump discontinuity, and T-C^{T") — 0. The set J"^ is iJ^-rectifiable 
(for the one-dimensional Hausdorff measure) and so, except for a set of H^- 
measure zero, can be covered by a countable family of graphs. In the 
sense, we can also define the approximate limits u±{t,x) from each side of the 
shock at {t,x), and A {t,x) denotes the shock speed. 

Lemma 3.1 (Characterization of the nature of shocks). Consider a scalar con- 
servation law with general flux / : R — s- K and a shock wave connecting the 
states U- , M+ and propagating at the Rankine-Hugoniot speed 



a{u-,u^) := / duf{su^ + {l — s)u+)ds. 



(3.2) 



For any real u' the nature of the shock wave in the averaged speed function 



a{t, x) = 



a_ := a(u_, u'), x<a{u^,u+)t, 
a.+ ■.= a{u+,u'), X > a{u-,u+)t, 



is uniquely determined by the signs of p(u^,u^,u') := (u_ — u') (u+ 
the jump (a-|_ — a_) of the averaged speed, as follows: 

1. If u± — u' have the same sign, 



(3.3) 



u_) and 



p{u-,u+,u) (a+ - a_) 



< 0, 5=^, 
> 0, T''. 



(3.4) 



2. If u± — v! have opposite signs, 



p(-u_,-u+,m') (a+ - a_) 



< 0, 

> 0, L"". 



(3.5) 



With some abuse of notation, in the first case of (|3.5p we have allowed here 
the case of equalities in the definition of a rarefaction-shock. 
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Proof. The inequalities can be checked geometrically from the graph of /, or 
else algebraically from the jump relation, as we now do. In view of 

a(u-, (w+ — M_) + a(u_, u') (m_ — u') — a(7i+, u') (u^ — u') 
= f{u+) - + - f{u'}) - {,f{u+) - f{u')) = 0, 

we obtain the identity 

(a{u-, u+) — a(u_, u')) (u_ — u ) ~ — (a(u_, u+) ~ a{u+, u')) (it+ — u'). (3.6) 

Let us recall Definition 12.11 which provides the classification of shocks. The 
identity (j3.6p shows that the shock in the function a is a compressive or rar- 
efaction shock if u± — u' have opposite signs; it is undercompressive if u± — u' 
have the same sign. (In fact, p.6p is nothing but the identity (|2.1ip already 
derived for general systems, but for scalar equations (|2.1ip reduces to saying 
that 7i+ 7r.) 

Next, it is interesting to observe that the identity p.6p can be cast in the 
new form 

(a(M+, m') — a(u_, u')) (u- — u ) = (a(u+, u') — a(w_, (u+ — w-). (3.7) 

This second identity allows us now to identify the nature of undercompressive 
shocks (either slow or fast) when v! — u± have the same sign, and the nature of 
compressive and rarefaction shocks when v! — u± have opposite signs. □ 

A result similar to Lemma [3.11 including the key identities (|3.6p and (|3.7p . 
will be established for fluid dynamics equations in the course of the proof of 
the following theorem and, later, a version valid for general systems will be 
established as well. 

Theorem 3.2 (Averaged matrix of scalar and fluid dynamics equations). Let u = 

u{t,x) and u' — u'{t,x) be solutions with bounded variation and arbitrary large 
amplitude to 

1. a scalar conservation law with general flux-function, or 

2. the Lagrangian p-system of isentropic fluid dynamics with general volume- 
pressure law p — p(y) satisfying solely the hyperbolicity condition dvp{v) < 
0. 

Then an entropy admissible shock of u or u' can not be a rarefaction- shock of 
the averaged matrix A.{t,x) :— A{u,u'){t, x) . 

Remark 3.3. 1. As observed earlier, the averaged matrix is not unique in 
general (except when N = 1). For the p-system we use the averaged matrix 
defined by the general formula (jl.5p . Later we will see (e.g. in Section\^ that 
for the Euler equations the averaged matrix must be based on the special structure 
of the equations. 

2. Theorem \3.2\ vrovides us with a key property required for L^ stability. Of 
course, approximate solutions u,u' (for instance constructed by front tracking) 
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generally admit non- admissible shocks; however, the strength of such shocks 
vanish in the limit (as the approximation parameters tend to zero) so that they 
do not prevent the stability estimate to hold in the limit. 

Proof. The statement for scalar conservation laws is consequence of Lemma |3. II 
and the entropy condition. Namely, consider a point {t, x) where the solution 
admits an entropy admissible shock (u-, u+) while u' is locally constant. Recall 
Oleinik's entropy inequalities |39| 

/(.)-/(z._)^/(^+) -/(._) for , between and u+. (3.8) 

V — U- u+ — U- 

According to Lemma 13.11 to exclude rarefaction-shocks we must exclude the 
first case in p.Sp where u' — u± have opposite signs and, simultaneously, 

{u -u-){a+-a-){u+-u-) <0. (3.9) 

Suppose (|3.9p hold with, for definiteness, w+ < u_. Hence w+ < u' < u^, so 
u+ — w_ < and u' — u- < 0. It follows from (|3.9p that a+ — a_ > 0. On the 
other hand, according to Oleinik's condition (|3.8p 

a > a+, 

which leads to a contradiction. Consequently, the first case in (|3.5p is impossible, 
and the shock is compressible or undercompressive, but can not be a rarefaction- 
shock. This completes the proof of the theorem for a scalar conservation law. 
Consider next the p-system 

dtu + dMv)^0, dtv-d^u = 0, (3.10) 

where u and v represent the fluid velocity and the (positive) specific volume, 
respectively. The pressure function p = p(v) is assumed to be decreasing, which 
ensures that p.lOp is strictly hyperbolic. Given two solutions {u,v), {u',v') of 
the p-system and setting ip u ~ u' , ip ■= v — v' , we obtain 

dtip - dx{c{v,v'f ip) = 0, dtLp-dxip = 0, 

where 

c{v,v ) = \ . 

V V ~ V 

Therefore, the averaged matrix (jl.Sp associated with the p-system is 

A{u,v-u',v') ^A{v,v') := (^_°^ -^C^J-^^')'^ ^ (3 

where \j(u,v]u,v') := ±c.{v,v') are the averaged wave speeds and correspond 
to the physical wave speeds \j{u,v) = -iic{v) := zLy^—dvp{v). 

Observe that the averaged speed of the p-system only depends upon a single 
variable, i.e. the specific volume v,v'. Considering a 2-shock wave connecting 
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(u_,w_) to {u+,v+) at the (positive) speed c(w_,w+) and denoting by {u',v') 
an arbitrary constant, we find the identity 



(3.12) 



(c(u_, v+) — c(w_, u')) {v' — w_) 
= k(u_, «+, u'){c{v-,v+) - c{v+, v')) {v' - v+), 

where 

k(?;_, «+,«') := {c{v-,v+) +c{v+,v')) (c(w_,'y+) +c(u_,'y')) ^ 
This shows that the 2-shock of the averaged matrix 

Ait,.) = \f-^'}; -<y-'-+);' (3.13) 

is a compressive or rarefaction shock if v± — v' have opposite signs, while it is 
undercompressive if v± — v' have the same sign. 
We can also derive the second identity 

(c(v+,v') — c(w_, v')) (w_ — v') 

, /- , - N (3-14) 

= k{v ,v+,v-) {c{v+,v ) " c{v^,v+)) {v+ - U-), 

which allows us to distinguish between compressive and rarefaction 2-shocks, 
on one hand, and, on the other hand, between slow or fast undercompressive 
2-shocks. This leads us to an analogue of Lemma [01 for the p-system: setting 



p('U_, «+, v') := (w_ - v') (w+ - w_), 

we have : 

1. \i v± — v' have the same sign, 



p(^_,^+,«')(c(«+,t'')-c(^;^,«')|jQ (3.15) 



2. If v± — v' have opposite signs, 



{< TZ^ 
> o' ^^'^^^ 

In turn, recalling that Wendroff's entropy criterion for 2-shocks of the p-system 

m 

c{v,V-) > c(w_,w+) for all v between v- and v+, (3.17) 

we can check that that the first case in (I3.16P can not arise if the shock is entropy 
admissible. We conclude that, for the p-system, admissible 2-shocks can not be 
rarefaction-shocks of the averaged matrix (|3.1ip . A similar conclusion holds for 
shock waves in the first characteristic family, and this completes the proof of 
the theorem. □ 
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3.2 Genuinely nonlinear systems 

In the case of genuinely nonlinear characteristic fields we can derive a mono- 
tonicity property of the eigenvalues of the averaged matrix. 

Lemma 3.4 (Shock speed monotonicity property). Let (jl.ip be a strictly hyper- 
bolic system of conservation laws with flux f : B^^ — s- admitting a genuinely 
nonlinear eigenvalue A^. Then there exists Si < 6q such that the following prop- 
erty holds for every averaged matrix satisfying p.4p . For arbitrary w_, u' g B^j^ 
the averaged wave speed Xi{-,u') is a strictly monotone function along the i- 
shock curve issuing from More precisely, if an i-shock {u^,u^) satisfies 

Lax shock inequalities 

Xi{u_) > Xi{u^,u+) > Xi{u+), (3.18) 

then the averaged speed satisfies the inequality 

X^{u^,u') >X,{u+,u'). (3.19) 

Furthermore, this holds for the p-system of Lagrangian fluid dynamics even for 
shocks with arbitrary large amplitude, provided the equation of state is genuinely 
nonlinear (that is, dy^p{v) > Oj. 

Proof. The right-hand state — u-|_(e) can be viewed as a function of u_ and 
a parameter e varying in the neighborhood of 0, with 

u+(e) = u- +er.,{u-) + O(e^). 

Then, we compute 

Xi{u+,u') - Xi{u-,u') = e ViAi(m_,u') • ri(u_) + O(e^) 

= e ViA,(w_, u_) • r,(u_) + 0{e^) + 0(e \u - u-|) > 0. 

Using now that Xi{u, u) ~ Xi(u) and the symmetric property^(u, u') —A{u', u), 
we obtain 2 ViA^ • = VA^ • > 0. Provided |e| + |u' — u_| < Si is sufficiently 
small, we conclude that the function Xi{-,u') is strictly monotone along the 
shock curve. 

If the shock (u_,m+) satisfies the entropy inequalities ()3.18|) and the nor- 
malization VAi • > is chosen, then e < and we conclude from the above 
calculation that the averaged speed decreases from u_ to w+. 

For the p-system, using the notation introduced earlier we can write 

dc 1 

dv c[v,v') 7o 

which keeps a constant sign if dyyp has a constant sign. Therefore, under the 
genuine nonlinearity assumption, the averaged speed of the p-system is strictly 
monotone along each Hugoniot curve. □ 
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We deduce that: 



Theorem 3.5 (Averaged matrix of genuinely nonlinear systems). Let u = 
u{t,x) and u' = u'{t,x) be (small amplitude) solutions with bounded variation to 
a system of conservation laws (jl.ip with genuinely nonlinear flux f , and let A be 
an averaged matrix satisfying (|1.4p . Then, an entropy admissible shock of u or 
u' can not be a rarefaction- shock of the averaged matrix A{t, x) :—A{u,u')(t,x). 

For scalar conservation laws and fluid dynamics equations, the monotonicity 
property exhibited in Lemma [3 . 41 provides us with a second proof of Theorem l3.2l 
which, however, requires genuine nonlinearity. 

Proof. We rely on Lemma 13.41 and consider a point (t, x) of jump of the BV 
function u: the limit U- is connected to u+ by an i-shock wave with speed 
Xi{u-,u+) satisfying Lax shock inequalities p.lSp . If the BV solution u' is 
continuous at (f, a;) then the inequalities characterizing a rarefaction-shock. 



3.3 Systems with general flux 

For characteristic flelds that are not genuinely nonlinear, a suitable version of 
Lemma 13.41 and Theorem 13.51 will now be established, by restricting attention 
to "robust" patterns corresponding to "strongly dominant" char- 

acteristic components, in a sense introduced below. 
We begin with a preliminary result. 

Lemma 3.6. Given an i-shock wave connecting the states G B^^, and 

given u' G B^q , define 



clearly contradict the monotonicity property (j3.19p . 



□ 




~lj{u±,u), 



A,(u_,'u+), r 




rj{u±,u') 




Then, the folllowing identities hold 




e 




(3.20) 



where the remainder Q, satisfies 



n^O{l)\r+ ^r,^ 



7l|A7-A||ar| + 0(e-)^|a-|. 
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In the course of the proof given below, we will derive the more general 
identities 



(A, - A, ) {u+ - u ) - (A - Aj ) (u+ - U-) 

= (A+ - AJ - u') + (A+ - A) {u+ - = 0(1) 

together with the following expression for the remainder 

r! : = ^ I (A - A+) a+ - (A - Xj) aj \ + O(e-) J2 \^J\. 



(3.21) 



(3.22) 



Proof. We can write 

(Ai(u+, u') - Ai(-u_, u')) (u- - u') + (Ai(m+, u') - Ai(w_, u+)) (w+ - w_) 
= (Aj(m+,u') - Xi{u-,u')) (m+ - u') - (Ai(u_,u+) - Xi{u-,u')) (m+ - u_) 

= ^, 

(3.23) 

where, using that u+ lies on the Hugoniot curve from u_ , 

h :=- (/(u+) - f{u') - X,{u+, u') {u+ - u')) + (/(ii_) - f{u') - K{u^,u') (u_ 
= - {A{u+,u') ~-J.i{u+,u')) {u+ - u') + (A(?i_,u') - Ai(ii_,M')) (ii_ - w') 
= - ^ (Ai(u+,u') - Aj (m+,m')) aj(u+, u') rj(u+, u') 

+ (Aj(w_, u') — Aj(u_, m')) aj(M_, u') rj{u-,u'). 
Therefore, introducing the notation 



"ff = 7^ u') := (A - Xf) af, 



the remainder takes the form 



Observing the decomposition 



a|-A,^ , , , , ^ /Xt -Xt . AT - A 



= 0(l)^|7+-77l+0(e^)El77L 
we arrive at the general statement p.2ip - p.22p . 
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Multiplying now (|3.23p by and if respectively, we obtain 

(X; - A) 1+ . (u+ - - (A+ - j;)a+ = -1+ ■ h, 
(Xi - A) ;r . [u+ - u_) - (A^ - A, ) = -ir . 

which gives p.20[) . Next, using u± — u' =: J2j ^^f find 

-(Af - A) ^ a+ r+ + (A - A^) ^ a+ r+ = h, 

3 3 

thus 

ifrf-liri = (A,^ - A) ^a+r+ + ( A - A, ) X!"7^7' 

which leads to the desired estimate p.20p along the same lines as in the proof 
of Lemma 12.51 □ 

The relations p.20p are the analogue, to the case of systems, of the identity 
(|3.7p found earlier for scalar equations. They allow us to relate the signs of 
Aj — A with the sign of A^ — Aj . This is so only when the reminder can be 
neglected, and this motivates us to propose the following definition. 

Definition 3.7. With the notation in Lemma \3.6\ the i- characteristic compo- 
nent across an i-shock is said to be strongly dominant if its left- and right-hand 
limits of satisfy 

|A+ - A7I |a±| > |r+ - rr| | A^ - A| |ar| + - ^ \al\. 

In this case, the shock is said to be robust relatively to the state u' . 

Following }26j we parametrize the i-wave curves with a global parameter 
fii{u) such that 

V/ii • ri{u) > 0. 

The hypersurfaces along which /i^ is constant determine a smooth foliation of 
the phase space and all i-curves are transverse to them. Observe that 

sgn(Ii(u_,M') • (u+ - = sgn(7i(M+,M') • (u+ - 

= sgn(//j(M+) - fii{u-)), 

when lies on the Hugoniot curve from u_ since this curve is tangent to the 
i-eigenvector at it_. 

From the identities in Lemma 13.61 we deduce immediately: 
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Lemma 3.8 (Characterization of the nature of robust shocks). With the nota- 
tion in Lemma \3.6\. the nature of a robust i-shock of the averaged matrix 



A{t,x) 



]A{u-,u'), X < A,;(m_, u+) 
IA(-u+,u'), X > Xi{u-,u+)t, 



is determined by the signs of pi{u^,u^,u') :— ai{u^,u') {fii{u-^-) — fii{u^)) and 
the jump (Ai(u_|_, u') — Ai(w_, u')) of the i-eigenvalue, as follows: 

1. If ai{u±,u') have the same sign, 

< 0, Sf^ and robust, 
Pi{u^,u+,u') {Xi{u+,u ) - Ai(-u_,u')) { (3.24) 



> 0, J-f^ and robust. 



2. If ai{u±,u') have opposite signs, 



> 0, Lf- and robust, 
p^{u_,u+,u') {Xi{u+,u') - Xi{u-,u'j) <( (3.25) 

< 0, TZf^ and robust. 

Remark 3.9. For scalar conservation laws every shock is robust (since D, = 0) 
and we can recover Lemma \3.1\ from Lemma[ 



Using the above result together with the entropy condition which for gen- 
eral system imposes that {u^,u+) is entropy admissible [311 [3S] if a-nd only if 
the shock speed Xi{u-, •) achieves at the point its minimal value along the 
Hugoniot curve, that is 

Xi{u-,u+) < Xi{u^,v) (3.26) 

for every state v along the Hugoniot curve between u_ to m+. Combining now 
the conclusions of Lemma 13.81 with the entropy inequality p.26p we arrive at: 

Theorem 3.10 (Averaged matrix for systems with general flux). Let u = u{t, x) 
and u' — u'{t,x) be (small amplitude) solutions with bounded variation to the 
strictly hyperbolic system of conservation laws (jl.ip with general flux f . Then, 
a robust and entropy admissible shock of u or u' can not be a rarefaction- shock 
of the averaged matrix A{t, x) :—A{u, u'){t, x). 



4 Stability for a class of linear hyperbolic sys- 
tems 

We are now in a position to derive an stability property for a large class of 
systems, by assuming the existence of a weight w = w(t, x) satisfying certain 
constraints that we specify. Consider a solution -0 — ''Pit, x) of a uniformly 
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hyperbolic system with bounded variation (|1.6p , together with its characteristic 
components aj , f3j defined earlier in Section [2l Relying on Lemma 12.31 on the 
evolution of the weighted norm, our goal is to establish that the right-hand side 
of (|2.7p is non-positive. At an i-shock. Lemma 12.51 and, more precisely, (|2.13p 
shows that the components and \(3'^\ coincide "up to first-order", with 
error terms proportional to | A+ — A^ | or j?-^ — rft_ \ ■ The weight w should have 
jumps that precisely compensate the effect of these error terms. In the present 
section, we assume that the matrix-field A satisfies the condition 



A 

1+ 



^^-l<|Af+-Af_| (4.1) 



at every jump point, and we formulate the jump conditions in term of the 
eigenvalues (A^ — A^) rather than the eigenvectors. 

Definition 4.1. With the above notation, the i- characteristic component across 
an i-shock is said to be strongly dominant if its left- and right-hand limits af 
satisfy 

«2 |A+-A7||a±|>|A+-A"|^KT|. 



We seek for conditions on — wJ ensuring that the right-hand side of (|2.7p 
is non-positive. In view of p.lSp . error terms associated with non-dominant 
components can be estimated by the dominant ones 

E 1/371 <G E \f^7l (4-2) 

j j dominant 

and no condition will be imposed on the jump — of non-dominant com- 
ponents. On the other hand, in view of the sign properties (|2.9p - (|2.10p of the 



More precisely, we impose that at each {t, x) G J^f^ and for j ^ i, 



characteristic flux Pj, it is natural to require that the weight w satisfies (I2.2ip . 



^ _ \ > K |A+ — A I, j < i and j dominant, 
< —K |A+ — A^l, j > i and j dominant 



w J — W^- 



(4.3) 



for a sufficiently large constant K > 0. 

Specifying the jump — of the i-component across an j-shock is more 
delicate. Observe that no condition is necessary on Cf- (compressive shocks) 
since both fiux have a favorable sign. On the other hand, it is hopeless to 
try to impose a condition on TZf- (rarefaction shocks) since both f3f^ have an 
unfavorable sign. Recall here Definition 13.71 and Lemma 13.81 which show that 
only strongly dominant rarefaction-shocks have been characterized, however. 
This motivates us to allow only such waves in the right-hand side of the key 
estimate (|4.6p below. 

Next, we focus on undercompressive shocks associated with dominant com- 
ponents. At this juncture, we recall Definition l3.7l and Lemma [3. 81 which provide 
us with the nature (slow or fast) of an undercompressive shock, but under the 
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assumption that the characteristic components are strongly dominant. This 
motivates the following constrains on the weight: 



_^ _ I < —K \Xf^ — Xfl], Sf^ and strongly dominant, 
* ' |> X|A^ — A^l, J^/^ and strongly dominant, 

while for dominant but not strongly dominant component we require the weaker 
condition 

I _ \ < K \\f\ — \f^\, Sf- and i dominant, , , 

w; -w, < ~ ' Y ' ' \ (4.5) 
\> -KlXf'^ - A^l, Tf" and i dominant. ^ ^ 

No condition is required on if the ^-characteristic component is not dominant. 

Proposition 4.2. Let A = A{t,x) be a uniformly hyperbolic matrix with 
bounded variation satisfying (j4.ip . Then, given any piecewise constant solution 
tp = il)(t,x) of p.6p and a piecewise constant weight w = w(t,x) satisfying the 
general conditions stated in Section\^ together with (|4.3p - (|4.5p . the associated 
weighted norm satisfies 

^llV'Wllw(t) < sup \^+{t,x)-^J-{t,x)l 

where ICf- C TZf- denote all i-rarefaction-shocks associated with strongly dom- 
inant characteristic components. Consequently, the solution ip satisfies the 
stability estimate 

U{t)\\LHR)<U{0)\\LHR)+ [ sup V"(s,a;)|ds, t e R+. 

Jo x£K^(s) 

(4.6) 

We conclude that the stability property (|1.7p holds when A does not 
admit strongly dominant rarefaction shocks. This result will be used in Section[6] 
in combination with our earlier conclusion (Theorem I3.10p that the averaged 
matrix of two entropy solutions of a system of conservation laws does not admit 
strongly dominant rarefaction-shocks. 

It is worth pointing out also that Proposition 14.21 remains valid for approx- 
imate solutions that satisfy the equations (jl.6p up to a measure source-term, 
provided the total mass of the source is added to the right-hand side of (|4.6p . 
See El] for details. 



Proof. We use Lemma 12.31 and estimate each term in the right-hand side of 
(EH): 

3 
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associated with a given i-shock. In the following, we will often make use of 
(|4.2p . Using the expansion (|2.13p we can write 



B=J2 {PJ w7 + sgn(/3+) 1/3- 1 



4) 

3 



+ 0{\rt, - rt \m I + 0(|A+ - A-|) J] |/5,7|, 

SO that by and dm]) 

S =(w,r sgn(/3-) + w+ sgn(/3+)) |/3r | _ ^(w+ - wj) \P; \ 

EK - ^t) \PJ\ + 0{\\t Af_|)|/5r| + 0(|A+ - A-|)^ 

Thanks to (|4.3p . for j ^ i the dominant j-components lead to a favorable 
sign of the jump — , while non-dominant components for j ^ i can be 
collected in the remainder using (|4.2p . We obtain 

B <(wr sgn(/3-) + w+ sgn(/3+)) |/?n 

- if |A+ - A-| 5] |/3-| + 0(|A|V - Af_|)|/3r| + 0(|A+ - A-|)^ |/3-|. 

The sum over dominant components allows us to suppress the remainder, pro- 
vided that K is sufficiently large (so that i4r|A+ — A^j dominates 0(|A+ — 
A~|)). We arrive at the inequality 

B <(wr sgn(/3-) + w+ sgn(/?+)) 1/3" | 

+ 0(|Af+ - AM) \f3-\ - I |A+ - A-| ^ 1/3,- 1 . (4.7) 

It remains to deal with the term \(3~\, which can be assumed to be dominant, 
since non-dominant components \(3~\ can be handled with ()4.2p . We will now 
distinguish between the cases of compressive, undercompressive, and rarefaction 
shocks. 

Case 1. If the i-shock is compressive, then we have sgn{P~) — sgn(/3+) = — 1 
and therefore 

B < -2u;-- |/3rl + 0(|At - AM) |/3-| - ^ |A+ - A-| ^ |/3-|. 



So we obtain 



S<_^min|^-|_||A+_A-| ^1/371 inrf. (4.8) 
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Case 2. If the z-shock is undercompressive, we will prove that 

S<-f |Af+-Af_||/3n - I |A+-A-| 5^1/371 inS^u:F^. (4.9) 

Here, we have sgn(/3~) = — sgn(/?j^). 

First of all, if the i-component is strongly dominant, then by (|4.4p we have 
determined the jump — in such a way that it can compensate the error 
term — |/3~|, provided the constant K is chosen to be sufficiently large 
so that K |A+ - A-| is larger than 0(|A+ - A"]). This leads us to (gl]) in the 
strongly dominant case at least. 

Second, when the z-component is dominant but not strongly dominant, we 
can assume (the other case being similar) that 

^2 \\t Af_| |a+| < \rt - rt\\f3-\ + |A+ - A-| ^ |/3,7|, 

thus by dUD and ([^1^ 

y \Xt - K\ < |A+ - A-| ^ (4.10) 

Here, the weight satisfies solely (14. 5p and we write 

(3- wr + (3+ w+ = sgn(/3-) - |/3+|) wr + sgn(/?-) (w" - w+) 
< K0{1) |Af+ - Xt\ 1/3-1 + 0(1 A+ - A-|) ^ 1/3, 



fe 



But implies 



KO{l) \Xt, Xt\W\ < I |A+ - A-| ^ 1/3,7 1. 

fc 

provided So , K2 are sufficiently small. In turn, taking into account the other 
characteristic components which yield a decay of X |A+ — A^| J2k \l^k\' 
arrive at 

B<~K\A+-A-\Y^ \f3-\ + I |A+ - A-| ^ 

k k 

= -||A+-A-|^|/37|, 

and we conclude that (|4.9p holds for all undercompressive shocks. 

Case 3. Finally, consider the case of an i-rarefaction-shock, for which no 
constraint has been imposed on the component w^. Here, we will show that 

B <C\A+ A-\ - - ||Af+ - \t\ |/3rl 

K ^ (4.11) 

- \A+-A-\J2\P~\ in7^f, 
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as well as the sharper estimate on non-dominant components: 

i? < - f I Af+ - AM 1/3- 1 - :| I A+ - A- 1 ^ 1/3- 1 in 7^f and non-dominant. 

(4.12) 

We first derive (|4.1ip for all rarefaction-shocks. From (|4.7I) we get 
B < 2w^^- |/3-| + 0(|At ~ \t\) W\ - y |A+ - A-| ^ |/3-| 

(4.13) 

< 3«;"- 1/3-1 --IA+-A-I 5] 1/371 

and we distinguish between two subcases: 
- If a,- a+ > 0, then by ((2?TT|) we have 

(A^-Af+)a+ + (Af_-A^)a- = 0(|Af+-AM)|/3-| + 0(|A+-A-|)^|/3-|. 

The two terms on the left-hand side above have the same sign, therefore 

1/3+ 1 + |/3n - I (A^ - At ) 1 + I (At - A^) a- \ 

= 0(|At - At|)|/3-| 0(|A+ - A-|) ^ |/3-|. 

We can then suppress the term 0(|A+ — A"|) above by taking K sufficiently 
large, and so 

i3<-||A+-A-| ^1/371, 
which -using once more the previous inequality- implies (|4.1ip . 

■A 



If Qij < 0, then using A^ < A < At wc find 

w\ = iAt - A^i i^n < I At - At I H - ".n 

<0(|A+-A-|)|7^+-V- 



(4.14) 



so that 



B < 0{\A+ A-|) - V'-l - f |A+ - A-| ^ 1/3^ 



where, in the latter, the estimate (|4.14p on |/3j | was used once more. Since 
the total variation X]2;GJ(t) l-^+(^' ^) ~ ^)l uniformly bounded in t, this 

proves (|4.1ip . 
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Finally, to treat the case of a rarefaction-shock that is not strongly dominant, 
we modify the argument in (|4.14p by using the condition (|4.10p : 



i/3n = iAf--A^iia,ri<iAt-Af-ii«ri 

<0(|A+-A-|)^|/3,-|, 



(4.15) 



which leads to 

The proof of Proposition 14.21 is completed. 



□ 



5 A new functional for nonlinear hyperbolic sys- 
tems 

5.1 Scalar conservation laws 

Conditions on the weight function w were proposed in Section |4l and our aim 
is now to show that such a weight exists. More precisely, we introduce now 
slightly different conditions which are easier to work with at this stage of the 
analysis and will be shown in the following section to imply the conditions 
required for our analysis of Section U) Importantly, we are going to see now 
that wave cancellations must be taken into account to establish that the weight 
remains uniformly bounded. The new technique proposed here relies on the 
adjoint system (|2.5p and on the property of generalized characteristics for the 
averaged matrix. The following is a natural generalization of the technique 
proposed earlier for genuinely nonlinear equations [251 131j which also rely on 
wave cancellation. 

For scalar equations, we use Dafermos' front tracking technique [13] . 

Theorem 5.1 (Functional for scalar conservation laws). Consider a scalar con- 
servation law with general flux f . Given Ci, K > there exists C2 = C2(Ci, K) 
such that the following property holds. Given (piecewise constant) front tracking 
solutions Uh,u'i^ satisfying the uniform bounds 



there exists a (piecewise constant) weight function w = w[u/j,uy satisfying 



together with the following conditions on undercompressive shocks of the aver- 
aged speed a := a{uh, u'jj: 



\\uh\ 



+ IKIk- +supTv^(uft(t)) + supry«(<)) < Ci, 



(5.1) 



I/C2 < w < C2 




5^ 



(5.2) 
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From this theorem, we deduce that the corresponding weighted functional 
is (essentially) decreasing: 

I \u'f^~ Uh\w[uh,u'f^]{t2)dx + I [A*^ - a~| |u_ - u'_\dt 

+ K |a+-a-||r-a-||M_-u'_|di (5.3) 

< / \u'y^- Uh\^[uh,u'^]{ti)dx + 0{h), t2>ti, 



where A denotes the shock speed of the function a. The remainder 0{h) is 
bounded by the maximum size of rarefaction fronts in Uh,u'f^, which tends to 
zero with the discretization parameter h. Hence, by letting /i — > in (|5.3p 
one deduces that the limits u := limw/i and u' := limu'^ satisfy the sharp 
contraction property 

/ \u' — u\w{t2) dx < / \u' — u\-w{ti) dx, t2>ti, 
Jr Jr 

where w := limw[M/i,My and the inequality is strict (except in the trivial case 
u = u'). Moreover, sharp dissipation terms can be made explicit, as was done 
for convex flux by Dafermos [13] and Goatin and LcFloch [20. . 

Proof. 1. For simplicity in the discussion and without loss of generality (modulo 
arbitrarily small perturbations of the data) , we can assume that all shocks in 
the function a satisfy the inequalities in Definition [2Tl] in a strict sense and that 
4> ■= Uh — u'f^ does not vanish, but of course may change sign at jump points. 
(This is possible since Uh,u'f-^ are piecewise constant with finitely many waves 
only.) 

Decompose the space-time M-|_ x M into finitely many maximal regions 
where Uh—u'y^ keep a constant sign and which, therefore, are limited by polygonal 
curves of changes of sign for the function t/i. According to Lemma 13.11 the 
boundaries of 17 consist of compressive or rarefaction shocks along which no 
constraint is imposed in (|5.2p . Consequently, we can turn attention to defining 
the weight w = w(t, x) within a given region fi. 

2. Consider backward generalized characteristics ym ■ [t^jtm] C [0, +oo) 
R associated with the averaged speed a and originating and finishing on the 
boundary of O, 

-£^{t) e I{a^{t,y„,,{t)),a+{t,y„i{t))) , t € 

where I{h, c) :— (min(6, c), max(6, c)). Observe that backward generalized char- 
acteristics within f2 are unique once their origin is given, since we assumed that 
no shock in a is left- nor right-characteristic. 



27 



More precisely, since u, v! admit finitely many waves and interaction points, 
we can select a complete family of generalized characteristics ym, rn = 1, . . . , M 
within the region 57, in the sense that every other characteristic meets precisely 
the same waves of Uh, ujj as one of the curves ym- 

3. We now determine the weight within the region by forma7iy solving the 
adjoint equation 

(9fW + adx'w = away from shocks, (5-4) 

in the following way. We prescribe the weight at the final point of each 
characteristic (when t„j < oo) 

w(Zm,?/,„(t„i)) = Co, m = l,...,M, 

where Co is a (large) positive constant. When < cx3 we require that w(i, y^n (t)) 
equals Cq for all sufficiently large times t. Then, we determine the value 
w(i, ym{t)) of the weight at arbitrary times t E (i.,„, im) by following the general- 
ized characteristic ym backward in time from the final time tm and by requiring 
that the weight: 

- remains constant away from shocks (so that (|5.4p holds), and 

- jumps up or down, according to the constraint (|5.2p . when it crosses a 
shock. 

Hence, our construction determines a weight defined along the curves ym 
and, in summary, we can write (away from shocks) 

w{t,ym{t)) = Co+K ^ 77,„(s) (a+(s,y„(s)) - a"(s,y„(s))), (5.5) 
se[t,tm) 

where the sign r]m{s) = ± is determined from (|5.2p and depend whether the 
shock is slow or fast undercompressive and whether it is crossed from left to 
right or from right to left. 

In turn, since the curves y^ is a complete family of characteristics, the weight 
can be uniquely extended as a piecewise constant function w = w(i, x) defined 
in the whole region SI. 

4. Next, for any fixed time T > we can introduce wave partitions (Liu [35j ) 
of the waves in both solutions Uh, u'j^ within the slab [0, T] x M. By definition, the 
solution Uh can be regarded as the superposition of finitely many elementary 
waves (wp,Mp-(-i) associated with piecewise affine functions ipp : [0,Tp] R. 
The latter determine the trajectories of the elementary waves in the plane. The 
wave (Mp,Up+i) is born at the initial time t = and completely cancelled out 
at the time t — Tp. The cancellation takes place at an interaction point with 
another wave. In the scalar case, no new wave is generated at interactions, and 
waves can only be split or cancelled out. A similar notation, {Ug,u'^^i) and 
ip'g : [0, T^] ^ R is used for the solution uj^. 
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5. We now establish that the function w remains uniformly bounded. It is 
(necessary and) sufficient to bound the maximum oscillation of w within the 
region 17 C [0, T] x M for arbitrary T and along each characteristic Um- We will 
prove that 

osc„i := sup_ •w{t,ym{t))- inf_ w{t,ym{t)) 

<rFK(o)) + rT/K(o)), 

where the implied constant depends only upon Ci, K arising in (|5.ip - ()5.2p . 

Recall first that fl contains only undercompressive shocks and that, within 
the region 17, the generalized characteristic ym crosses undercompressive shocks 
transversally. To derive (j5.6p we need to take into account wave cancellations 
and to rely on the following linearity property of the averaged speed a with 
respect to elementary waves: if a wave u+) in the solution u/,, is split into 
elementary waves (up,Mp+i) then for every constant u': 

a{u+,u') — a(u_ , u') — a{up+i,u') — a{up, u'), (5-7) 
p 

which we refer to as the linearity property. (Of course, this property holds for 
arbitrary functions.) 

To estimate the weight w along a path ym we fix an interval [ti, t2\ C (t_m,tm) 
and consider the set of all elementary waves (up, Up+i), p G E = Em{ti,t2), and 
{u'^,u'^^j^), q <E F = Fm{ti,t2), that are crossed by the characteristic ym at 
least once (but possibly more than once). Observe that a slow shock can only 
be crossed from right to left as we move backward, while the opposite is true 
for a fast shock. At this juncture we need to recall Lemma 13.11 which implies 
that a specific sign rjp = ± can be associated with a given elementary wave 
(up,Mp+i) which only depend whether Uh — u'f^ is positive or negative in the 
region fi, whether the wave is increasing or decreasing. (Note in passing that 
signs are different also for the waves (up, itp+i) and the waves (Up, Wp+i).) Then, 
at a point of discontinuity (s,j/m(s)) the sign rim{s) arising in (|5.5p is ±7yp and 
alternatively positive and negative as the wave is crossed from left to right and 
then move back to the left-hand of the characteristic at a later time. 

Consequently, given a pair (p, q) G E x F, together with the trajectories 
iPp,ip'q the contribution of the waves (up,Up+i) and (Wq,u^_).]^) to the weight 
along the curve ym can be estimated by the strengths of these waves 

\Up+l - Up\ + \u'q+i - u'g\ 

plus the alternating sum 

ihM) 

:= ^±(^(a(up+i,M^+i) - a('Up,u^+i)) - (a(up+i, u'^) - a(Mp, u^))) , 

over all pairs of waves crossings with the characteristic y„i within the time 
interval [ti, t2\. The sign is alternately positive and negative as one moves along 
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the characteristic. Consequently, we can keep only one term at most 



< 



|Up+l - Up\ - u'^\. 



In turn, by taking into account all pairs (p, q) of elementary waves we con- 
clude that 

osc™ < ^ \Up+l - Upl + ^ - Uql + X! '""P+l ~ ~ ""9! 

< TViu^m + rFK(o)) + rFK(o)) rFK(o)), 

which under the assumptions (jS.ip implies (|5.6p . □ 



5.2 Nonlinear hyperbolic systems 

We now turn attention to general systems (jl.ip and to piecewise constant so- 
lutions Uh generated by wave front tracking. For the actual construction and 
properties of front tracking approximatiorlj we refer to the general theory in 
[3, .251 [3H] (for general systems via the Glimm scheme) as well as to the ear- 
lier references [131 (THl SOI [SI O [ST] (for genuinely nonlinear systems via front- 
tracking) . 

Specifically, given two approximate solutions Uh,u'f^ we introduce the aver- 
aged matrix A ^ A(uh,u'i^) and the function Tp — u'/^ — Uh- Then, by (|2.3p we 
define the characteristic components aj and their traces = aj{t^x±) at a 
shock {t,x). Without loss of generality we can assume that shock trajectories 
in Uh,u'f^ cross at finitely many points only and never superimpose on an open 
interval of time, and that the components aj never vanish while the shocks 
are non-characteristic on both sides. This can be achieved by arbitrarily small 
perturbations of the data. Following [SS] we define the strength e^^ = e"^ {t, x) 
at a shock of the matrix-valued map A as the strength of the corresponding 
shock in the solutions u/j, u^. The latter is measured with the global parameter 
introduced in Section [5T51 

^^A^^ ^\^^^{ut) " ^M{u'^)\, i-shock in u,,, 

- i-shock in m'^. 

Clearly, we have 

|A+-A"|<e'^ at shocks. 

We now introduce conditions on the weight w = w(i, x) which are most 
convenient to work with at this stage. Recall first that the relevance of the 
signs of the characteristic components has been emphasized in (|2.19p - (|2.20p . 
and in Lemma 13.81 In particular, the latter allows us to decompose the set of 

^See the recent preprint: O. Glass and P.G. LeFloch, Nonlinear hyperbolic systems: Non- 
degenerate flux, inner speed variation, and graph solutions, Arch. Rational Mech. Anal., to 
appear. 
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all j-discontinuities in into two sets according to whether the quantity 

(defined at a jump point {t, x) from the left- and right-hand traces) 



Pi(w,, ,u^,u^), shock in M^, 



PiWh ^'u'h^,Uh), shock in 



is negative or positive respectively. 

J^^ = {p^{t,x) < 0} U {p^{t,x) > 0}. 

Lemma [3.8l shows some compatibility between the signs of the jump of the weight 
w^ — w~ and the jump of the eigenvalue Xi{u'^ ,u'f^'^) — Xi{uJ^ ,u'f^^) = A^ — A^. 

Actually, the above is true for strongly dominant shocks only. We pro- 
pose here to relax this condition and to impose the conditions suggested by 
Lemma 13.81 even for non-strongiy dominant components. We will see that this 
strategy leads to a well-defined and well-behaved weight w. 

In other words, we require that for all j ^ i 

w+-w- = / J<^anda7a+>0, 
^ ^ [-Ke^, j > i a.nd aj a+ > 0, 

while for j — i and all undcrcompressive i-shocks 



-X(A^-A^), Pi < and a7a+ > 0, 
K (A^ - A^), Pi > and af > 0, 



(5.9) 



which corresponds to well-defined conditions for each increasing/decreasing shocks 
in regions where ai is positive/negative. The proof of the following theorem will 
rely on a wave partition into increasing/decreasing elementary waves. 

Theorem 5.2 (L^ functional for systems of conservation laws). Consider a 
strictly hyperbolic system of conservation laws with general flux f . Given Ci, K > 
there exists C2 = C2{Ci, K) such that the following property holds. Given 
(piecewise constant) front tracking solutions Uh, wjj satisfying the uniform bounds 

\\uh\\L^ + \\u'Jl^ + supTV{uh{t)) + snpTV{u',M < Cu 
t t 

there exists a (piecewise constant) weight function w = (wj [li/j, u'^])i<j<7v sat- 
isfying 

together with the conditions ()5.8p and ()5.9p . 

Proof. 1. Our construction will be similar to the one proposed for scalar equa- 
tions, with however some important modifications due to 

- the new waves generated at interactions, and 

- the nature of shocks (fast/slow, compressive/rarefaction) which is not en- 
tirely determined by the signs of the characteristic components af and the 
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coefficient pt. Indeed, Lemma is concerned only with dominant components 
aj {j ^ i) and strongly dominant components Ui across i-shocks. 
However, by taking into account contributions due to small waves the construc- 
tion can still be carried out and wave cancellations be taken into account, as 
we now explain. Fixing some index j we focus on defining the component Wj 
of the weight. 

2. Decompose M+ x M into maximal regions Vtj where the characteristic 
component Uj keeps a constant sign and which are limited by polygonal curves 
across which aj changes sign. According to Lemma [3.8l the boundary of a region 
Vij consists of 

- fc-shocks {k 7^ j) with dominant characteristic components (but oppo- 
site signs), or 

- compressive or rarefaction j-shocks with dominant characteristic compo- 
nents (with opposite signs), or 

- shocks with non-dominant components (with opposite signs), 

- or possibly an interval of the initial line t — Q. 

We observe that along of the first type of boundary the simpler constraint (|5.8p 
is imposed while along the other types of boundary no constraint is imposed on 
the j-component of the weight Wj. 

Consequently, we need not only to deal with Wj in one of the regions fij 
(in which the constraints (|5.8p - (|5.9p are relevant), but also to take into account 
waves leaving one of the regions fij. 

Moreover, in contrast with the case of scalar equations, a region U,j may now 
contain not only undercompressive j-shocks but also 

- /c-shocks (fc ^ j) and 

- compressive shocks or rarefaction-shocks of the j-family (which however 
must be non-dominant). 

3. We will rely on a sufficiently large family backward generalized j-charac- 
teristics : [0,cx)) R associated with the eigenvalue 

^(t) e l{Xf_it,ylit)),Xf^it,yi^it))), t e (0,00). 

Here, a given region flj may contain compressive j-shocks so that backward 
generalized characteristics need not be unique. Also, such regions may well 
allow j-waves to cross their boundary: we impose that j-characteristics never 
exist a region f2j, except along a boundary consisting of a /c-wave with k ^ j. 
Therefore may only exit to the right of Vlj by crossing a fc-wave with k < j, 
or exit to the left by crossing a fc-wave with k > j. And a j-characteristic does 
not cross compressive and rarefaction j-shocks, but may travel with the same 
location and speed as a rarefaction shock on some time interval. 

Since the solutions Uh,u'i^ admit finitely many waves and interaction points, 
we can select a complete family of generalized j-characteristics , m — 1 , . . . , M 
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covering the whole of IR+ x M, such that every other j-characteristic meets pre- 
cisely the same waves of u,u' as one of the curves in the family y^. 

4. Next, backward from the time T > we can determine wave parti- 
tions within the slab [0,T] x M of the waves in both solutions Uh,u'^ (Liu [55]). 
By definition of a wave partition, the solution Uh is regarded as the linear 
superposition of finitely many elementary k-waves {Up,Up^i) with trajectories 

■ [Tp,Tp] —f M. The wave iUp,Up_^_l) is born at the initial time and 

is completely cancelled at the time t = T^. For systems, new waves may be 

generated at wave interaction points. The total strength, and the change along 

, 

their trajectories, of completely cancelled waves (Tp < T), new waves (0 < Tp), 


and surviving waves Tp — T) can be estimated uniformly by the initial total 
variation, the interaction potential, and the cancellation measure, all of them 
being in turn controled by the total variation of the initial data u/i(0). A sim- 
ilar notation, {u' q,u' ^j^i) and (^'^ : [0,T/] — > R will be used for the elementary 
fc- waves in the solution uj^. 

5. We now determine the weight by formally solving the adjoint equation 

9tw + A dx'w = away from shocks, (5.10) 

in the following way. We prescribe the weight at some sufficiently large time T 
along each characteristic 

w,(T,y^(T))-Co, m=l,...,M, 

where Co is a (large) positive constant. Then, we determine the value Wj(t, y4(i)) 
of the weight at arbitrary times t < T by following the j-characteristic y-^ back- 
ward in time and requiring that Wj 

- remains constant away from shocks (so that (|5.10p holds), 

- jumps up or down, according to the constraint (|5.8p when it crosses a 
/c-shock (fc ^ j), and 

- jumps up or down, according to the constraint (|5.9p when it crosses an 
undercompressive j-shock. 

By construction compressive and rarefaction j-shocks are never crossed by j- 
characteristics. 

Our construction determines a function defined along each of the curves , 
and away from shocks we can write 

ag[t,T) 

+ K 77„(s)e^(s,y^(s)), 

se(t,T) 

fc-shocka 

where the sign r]m{s) = ± is determined from (|5.8p - (|5.9p and depends whether 
the shock is slow or fast undercompressive and whether it is crossed from left 
to right or from right to left. 
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In turn, by using that the curves yl^ form a complete family of characteristics, 
the weight is uniquely extendable as a piecewise constant function — v^jit, x) 
defined in R+ x R. 

6. We now claim that the function Wj remains uniformly bounded, indepen- 
dently of the parameter h. We need to estimate the maximum oscillation of Wj 
in [0,T] X R along each ^'-characteristic t/^. We will actually prove that 

osc^j := sup J {t,y^^{t)) - inf Wj (t,y^(t)) 
te[o,T] *e[o,T] 

<TV{uhm+TV{u'^m. 

We follow here the argument from the scalar case, but now the regions flj may 
contain not only undercompressive shocks (which the generalized characteristic 
2/„i crosses transversally) but also fc-shocks and other (non-dominant) j-waves. 

The linearity property with respect to elementary waves, pointed out for 
scalar equations, is now expressed in terms of the j-eigenfunction Xj = ^f- if 
a j-wave {u_,u^) in the solution is split into elementary waves (M^,Up_|_]^) 
then for every constant u': 

Xj{u+, u) — \j{u^ , u) — Aj(u^^]^, u) — ^j{u-'p, u). (5.13) 

p 

To estimate the weight Wj along a given path y-j^ we consider an arbitrary 
interval [^1,^2] C [0, T] and consider the set of all elementary waves crossed by 
the characteristic y^^ at least once (but possibly more than once), precisely: 

- all elementary j-waves (w^,Up_|_i), P € E^^ = El^{ti,t2), associated with 
the solution Uh, 

- all elementary fc-waves {Up, Up_^_i), p e E^^ — E'^{ti,t2), with k ^ j, 

- all elementary j-waves (u^"* , u'^^j^), q £ F^-^ — F,n{ti,t2), associated with 
the solution u'^, and 

- all elementary fc-waves {u'^^'^ , u'^^^) , q G ~ F„i{ti,t2), associated with 
the solution w'^. 

Observe that a slow j-shock can only be crossed from right to left as we 
move backward along a j-characteristic, while we have the opposite direction 
for a fast shock. Relying on the key observation in Lemma we observe that 
a definite sign 77^ = ± (depending whether the wave is increasing or decreasing) 
can be associated with each elementary wave (m^, Up^-^), so that the sign 77m(s) 
in (|5.1ip is ztrjp where the sign ± only depends 

- whether the characteristic at the time s lies in a region Qj where the 
characteristic component is positive and —rjp in a region where it is negative, 

- whether the wave is crossed from left to right, or from right to left. 
Consequently, as the characteristic crosses through the elementary waves 

of Mfi , , alternatively positive and negative signs arise in the first sum in (|5.1ip . 
Indeed, this is true in a given region flj, as well as when the wave crosses its 
boundary and enter another region of constant sign. 



34 



Consequently, given a pair {p,q) € E^^ x F/„, together with the trajectories 
<p^, ip'-'q the contribution of the j-waves (u^, Up+i) and (u';^, u'^+i) to the weight 
Wj along the curve within the time interval (ii, is estimated by 

Iwp+i-^pl + l^'^+i""'?! 

plus the sum 

over all pairs of crossings of the characteristic t/m within the time interval \ti , ■ 
The signs alternate as a given wave passes from left to right, or then back from 
the right to the left of the characteristic, and we conclude that terms cancel out 
two at a time and a single term only must be kept at most 



< 



< -u^l -u\j 



In turn, by taking into account all pairs (p, q) of elementary waves and by 
making use of properties of wave partitions (uniform bounds on cancellation 
and interaction measures) we find 



' 9+1 13 1 

+ Y l"p+i-"plK+i-<l 

< TV{u,M) + TViu'^m + TV{uum TV{u',M)- 

On the other hand, to handle the second sum in (|5.11[) we observe that the 
^-characteristic ym crosses a fc-shock with k ^ j at most once, and have suitable 
alternating properties. This second sum accounts for the contribution 



<TV{unm + TV{u'^{Q)). 
The time interval (ii, being arbitrary we conclude that 

osc^„ < Tt/K(0)) + TFK(O)) + TViunm T1/K(0)), 
which shows (|5.12p and completes the proof of Theorem 15.21 □ 
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6 Continuous dependence property 



We are now in a position to prove the estimate (ll.2p . For definiteness we 
state our results for solutions constructed by front tracking (see Section [5^ . 
With some modification, our method also applies to solutions constructed by 
the Glimm scheme. 

Theorem 6.1 {L^ continuous dependence of solutions). Consider a strictly 
hyperbolic system (jl.ip with general flux f , under the assumption that there 
exists an averaged matrix A satisfying (|1.4p and 

\r,{u+,u')-r,{u^,u')\ < \K{u+,u') ~K{u-,u% (6.1) 

for all states u± , u' under consideration. Then, wave front tracking approxima- 
tions Uhyu'f^ with sufficiently small amplitude and total variation satisfy, for all 
times t > 0, 

IK(t) - uh{t)\\mR) < IK(o) - uh{o)\\mR) + o{h) (6.2) 

and, consequently, the limit solutions u = lini/j^o Uh and u' := lim/j^o ""Ji satisfy 
the continuous dependence property 

Wu'it) - u(t)|Ui(R) < \\u'{0) - u(0)|Ui(K). (6.3) 

This is a generalization of a theorem established earlier [25l [31] for the class 
of genuinely nonlinear systems. Indeed, for such systems, the condition (|6.ip 
always holds since, on one hand, 

\ri{u+,u') -ri{u^,u')\ < \u+ - 

and, on the other hand by Lemma 13.41 

|m+ — u_| < \Xi{u^,u') — Aj(7i_, u')|, genuinely nonlinear fields. 

We conclude with: 

Corollary 6.2. The continuous dependence properties (|6.2p (approximate so- 
lutions) and ()6.3p (exact solutions) hold for 

1. the Lagrangian p- system of fluid dynamics, and 

2. the Euler equations of isentropic fluid dynamics, 
with general pressure-density equations of state. 

Proof of Theorem \6.1l We rely on the framework developed in previous sections. 
Let w = (wj{t,x))i<j<N be the weight determined in Theorem 15.21 We need 
to check that the conditions (|5.8p - (|5.9p imply the conditions (|4.4p - (|4.5p required 
for the stability theory of Section [H Once this is checked, the desired es- 
timate ()6.2p follows from Proposition 14.21 (L^ stability for linear systems) and 
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Theorein l3.10l (non-existence of strongly dominant rarefaction-shocks), and the 
properties of front tracking approximations, especially the fact that rarefaction- 
shocks in Uh, u'l^ have maximal strength converging to zero with h. 

To compare the conditions in Sections 2] and [51 we use the classification in 
Definition 12.11 and the characterization (obtained in Lemma 13. 8p of the shocks 
in terms of: 

- the sign of characteristic components aj and 

- the sign of the jump of the eigenvalues of A —A{uh, u'fj. 
First of all, in view of (|2.19p {aja'j > for dominant j-components of an 



i-shock) the conditions (|5.8p imply 

Second, concerning strongly dominant i-components , we see that (|5.9p in 
combination with (|3.24[) show that 

Sf- and strongly dominant, 
Tf^ and strongly dominant, 

which implies ()4.4p . 

Third, for an undercompressive z-shock whose i-component is dominant not 
strongly dominant, the sign of the jump — is not well defined, and the 
sign of the jump — may not agree with the general rule (j2.2ip . However, 
from (15.911 we can still deduce that 




|w+ - wr I ^ iiT I A^ - A^ 1 , 5f U J^f- and i dominant, (6.5) 
which implies ()4.5|) . □ 



Proof of Corollary \6.2\ From (|3.1ip we see that the eigenvectors of the averaged 
matrix associated with the p-system 



- f i\ ( ±c{v,v') 



only depend upon the eigenvalues ±c(w, v'). This shows that (|6.ip holds for the 
p-system. 

Next, we turn attention to the Euler equations of isentropic fluid dynamics, 

dtp + dM^O, dtipu) + d^pu^ + p{p)) = 0, (6.6) 

where p > and u represents the specific density and velocity of the fluid, 
respectively. The pressure p = p{p) is assumed to be strictly increasing for 
p > 0, so that the equations are strictly hyperbolic, away from the vacuum 
{p = 0) at least. We use here the density p and the momentum q :— pu, 
which are the conservative variables of (j6.6p . To define the averaged matrix we 
consider two solutions {p,q), {p',q') of (|6.6p and write 



dtip-p') + d,iq- q') = 0, dtiq - q') + dj^ - ^ + p{p) ' pip)) = 0, 
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which we now transform into the form p.6p . 

Let us introduce a function e = e(/9, g, p', g') satisfying 



9^ _2 



that is, 



P P' 

e{p,q,p',q') 



e\p~p') + 2e{q-q'), 
pu — p'u' — \J pp' [u — u') 

p- p' 



(6.7) 



2 2^4 
Introduce also the function c : ^ M by 

c[P,P)-= ; • 6.8) 

P- P 

Then, the maps := p — p' and :— q — q' satisfy the hnear hyperbohc 
system 

dt^ + d,^ = 0, 9tV + 9,(-eV + 2e^P + c\p, p') ^) = 0, 
which corresponds to the following choice of averaged matrix: 

,2 4). (6.9) 

Clearly, the eigenvalues of the matrix A{p, g; p'; q') are Xj = e ± c, while the 
eigenvectors are Vj = (l,e±c)-^. This establishes the desired property (|6.ip for 
the Euler equations. 

Consequently, Theorem 16 . 1 1 applies to both the Lagrangian and the Eulerian 
formulations of the fluid dynamics equations, and this establishes (|6.2p and (|6.3p 
for the fluid dynamics systems. □ 

7 Concluding remark 

The construction of the weight proposed in Section [5] can be generalized, by 
noticing that the jump of the weight need not be related to the averaged speed 
a. 

Proposition 7.1. Let tt : ^ M 6e a given smooth junction. Under the 
assumptions of Theorem \5.1\ there exists a weight w satisfying, instead of (|5.2p . 
the more general condition 

^+ ^_ _ i-K{TT+ -TT-), -^p-)'^{J+ >0 and^p~'ilj+ >0. 



where ip := Ufi — uj^ and tt^ — n{u'^,Uj_'). 
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The proof is immediate from the one of Theorem 1 5. II Of course, we recover 
the previous setting by choosing tt = a. 

As far as the application to hyperbohc conservation laws is concerned and in 
order to use this weight in the weighted norm however, we need (|2.2ip to hold 
at undercompressive waves, that is in the region where i/' keeps a constant sign. 
This essentially forces us to pick up tt = a in Proposition 17. II 
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